Fuente: artículo online
library(sf)
## Linking to GEOS 3.9.3, GDAL 3.5.2, PROJ 8.2.1; sf_use_s2() is TRUE
library(mapview)
## The legacy packages maptools, rgdal, and rgeos, underpinning the sp package,
## which was just loaded, will retire in October 2023.
## Please refer to R-spatial evolution reports for details, especially
## https://r-spatial.org/r/2023/05/15/evolution4.html.
## It may be desirable to make the sf package available;
## package maintainers should consider adding sf to Suggests:.
## The sp package is now running under evolution status 2
## (status 2 uses the sf package in place of rgdal)
baseurl <- "http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?"
wfs_request <- "request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json"
fi_regions_wfs <- paste0(baseurl,wfs_request)
fi_regions <- st_read(fi_regions_wfs)
## Reading layer `OGRGeoJSON' from data source
## `http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json'
## using driver `GeoJSON'
## Simple feature collection with 121 features and 4 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -66.17672 ymin: -28.02735 xmax: -64.48182 ymax: -26.05595
## Geodetic CRS: WGS 84
mapview(fi_regions, label = fi_regions$name, color = "darkgreen", col.regions = "green", alpha.regions = .05)
Con gusto te explico el código!
library(sf)
library(mapview)
library(sf): Esta línea carga la librería sf, que es fundamental para trabajar con datos espaciales en R. sf nos permite leer, manipular y analizar datos geográficos de manera eficiente.
library(mapview): Cargamos la librería mapview para crear mapas interactivos de manera sencilla. Con esta librería, podremos visualizar nuestros datos espaciales en un mapa web. Construyendo la URL del servicio WFS
baseurl <- "http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?"
wfs_request <- "request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json"
fi_regions_wfs <- paste0(baseurl,wfs_request)
baseurl: Aquí se define la dirección base del servidor WFS. Un servidor WFS (Web Feature Service) es un servicio web que permite acceder a datos geográficos en formato vectorial.
wfs_request: Esta parte especifica la solicitud que se le hará al servidor WFS. En este caso, estamos solicitando un conjunto de características separados por & request=GetFeature: Indica que se solicita obtener datos de una capa (feature) del servicio WFS. service=WFS: Especifica que se trata de un servicio WFS. version=2.0.0: Define la versión del servicio WFS que se está utilizando (en este caso, la versión 2.0.0). typeName=espacio01:tuc_juris: Define el nombre completo de la capa que se desea descargar. (“espacio01:tuc_juris”) outputFormat=application/json: Indica que se desea recibir los datos en formato JSON.
fi_regions_wfs: Finalmente, se combina la dirección base con la solicitud para formar la URL completa del servicio WFS.
fi_regions <- st_read(fi_regions_wfs)
## Reading layer `OGRGeoJSON' from data source
## `http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json'
## using driver `GeoJSON'
## Simple feature collection with 121 features and 4 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -66.17672 ymin: -28.02735 xmax: -64.48182 ymax: -26.05595
## Geodetic CRS: WGS 84
st_read(fi_regions_wfs): Esta función de la librería sf se utiliza para leer los datos espaciales desde la URL del servicio WFS. Los datos se leen y se almacenan en un objeto de tipo sf, que es la forma estándar en R para representar datos espaciales.
mapview(fi_regions, label = fi_regions$name, color = "darkgreen", col.regions = "green", alpha.regions = .05)
mapview(fi_regions): Esta función crea un mapa interactivo utilizando los datos espaciales almacenados en fi_regions. label = fi_regions$name: Cada polígono en el mapa se etiquetará con el valor correspondiente de la columna “name” en los datos. color = “darkgreen”, col.regions = “green”, alpha.regions = .05: Estos argumentos se utilizan para personalizar la apariencia del mapa, definiendo el color del borde de los polígonos, el color de relleno y la transparencia del relleno, respectivamente.
Este código completo realiza las siguientes tareas:
Conecta a un servidor WFS: Se establece una conexión a un servicio web que proporciona datos geográficos en formato vectorial.
Descarga datos: Se descargan los datos de una capa específica (en este caso, “espacio01:tuc_juris”) del servidor WFS.
Crea un objeto espacial: Los datos descargados se convierten en un objeto espacial en R, lo que permite realizar análisis espaciales.
Visualiza los datos: Los datos espaciales se muestran en un mapa interactivo, donde cada polígono está etiquetado con un nombre y tiene un estilo de color específico.
En pocas palabras: este código te permite visualizar datos geográficos de un servidor remoto en un mapa interactivo utilizando R.